Quantum Andreev effect in 2D HgTe/CdTe quantum well-superconductor systems 
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The Andreev reflection (AR) in 2D HgTe/CdTe quantum well-superconductor hybrid systems is 
studied. A quantized AR with AR coefficient equal to one is predicted, which is due to the multi- 
Andreev reflection near the interface of the hybrid system. Importantly, this quantized AR is not 
only universal, i.e., independent of any system parameters and quality of the coupling of the hybrid 
system, it is also robust against disorder as well. As a result of this quantum Andreev effect, the 
conductance exhibits a quantized plateau when the external bias is less the superconductor gap. 

PACS numbers: 74.45.+C, 85.75.-d, 73.23.-b 



Recently, the topological insulator (TI), a new state 
of matter, has attracted a lot of theoretical and exper- 
imental attentionji"— The TI has an insulating energy 
gap in the bulk states, but it has exotic gapless metallic 
states on its edges or surfaces. The TI is first predicted 
in two-dimensional (2D) systems, e.g. the graphene and 
HgTe/CdTe quantum weU (QW)^i^ The 2D TI has the 
gapless helical edge states and exhibits the quantum spin 
Hall effect. This helical edge states, with the opposite 
spins on a given edge or opposite edges for a given spin 
direction containing opposite propagation directions, are 
topologically protected and are robust against all time- 
reversal-invariant impurities. Soon after that, people also 
found the TI in three-dimensional (3D) materials, e.g., 
Bii-xSbx, Bi2Se3, etci^ On the experimental side, the 
TI in HgTe/CdTe QW, Bii^^Sb^, Bi2Sez, and so on, 
have been successfully confirmedj^i^ii 

The Andreev reflection (AR) which was found about 
fifty years ago is an important transport process;^ The 
AR occurs near the interface of a conductor and a super- 
conductor, in which an incident electron from the metal- 
lic side is reflected as a hole and a Cooper pair is created 
in the superconductor. When the bias is smaller than the 
superconductor gap, the conductance of the conductor- 
superconductor hybrid device is mainly determined by 
the ARjiS For a perfect conductor-superconductor inter- 
face, the probability of the AR can reach one. However, 
the AR coefficient is usually very small due to various 
scattering mechanisms such as the contact potential of 
the interface, the impurities, the mismatch of the density 
of states of the conductor and superconductor, and so 
on. Although the AR has been extensively investigated 
in various conductor-superconductor hybrid systems, up 
to now, people has not observed a robust quantized AR 
that holds for various system parameters and sustains in 
the presence of a variety of scattering mechanisms. 

In this letter, we study the AR in a 2D TI- 
superconductor (TI-S) hybrid system. We found that 
in this system the AR is quantized with AR coefficient 
being one, because of the multi- Andreev refiection along 



the TI-S interface. Importantly, this quantized AR is 
universal and robust. The quantization of AR persists 
regardless of the system parameters such as the Fermi 
energy, energy of the incident electron and the size of the 
system. It is also robust against the presence of impu- 
rities and contact potential of the interface. Due to the 
quantum Andreev effect, the conductance exhibits the 
plateau with the plateau value being Ae^/h (2e^//i) for 
the two-terminal (four-terminal) TI-S hybrid system. 

Since the TI phase in HgTe/CdTe QW has been exper- 
imentally realized j^ii we shall focus on the HgTe/CdTe 
QW in the following calculation. The results will 
be the same for other 2D TI systems. We consider 
two HgTe/CdTe QW-superconductor hybrid devices as 
shown in Fig. 1(f): one is a ribbon of HgTe/CdTe QW 
coupled to a semi-infinite superconducting lead which is 
referred as the two-terminal device and the other is a 
cross of HgTe/CdTe QW coupled to a superconducting 
lead which is referred as the four-terminal device. The 
advantage of the four-terminal device is that the "trajec- 
tories" of the incident electron and reflected hole can be 
clearly shown. 

The hybrid devices are described by the Hamiltonian 
H — Hti + Hs + He, where Hti, Hs, and He are 
the Hamiltonian of the HgTe/CdTe QW, superconduct- 
ing lead, and the coupling between them, respectively. 
By discretizing spatial coordinates of the continuous ef- 
fective Hamiltonian Hti in Ref|3 and using the Nambu 
representation, the Hamiltonian Hti is given byJ^ 

i/„ =^*t^ii^rj+ ^ ¥,H;s^i+3 + H.C., (1) 

* a — {Sx,5y),i 

where i — {ix,iy) is the site index, Sx and Sy 
are unit vectors along x and y directions, — 
(oi, &i, C; , dj )"^, and ai,bi,Ci,di are annihilation oper- 
ators of electron on the site i at the states |s,t), 
\Px + iPyA), \s,i), and -|p^ - ipy,i). In Eq.(l), 

Hn/iSx/iSy - ( q''"''^"'' \ . - ) are the 4 X 4 
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tgp ~ A/2a. Here Ep is the Fermi energy (pinned by su- 
perconductor condensate), a is the lattice constant, and 
A, B, C, D, and M are the system's parameters which 
can be experimentally controlled. The Hamiltonian Hg 
of the superconducting lead is: Hs = Sk.o-ekOgkcr'^Skcr + 
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Here 



EkA(agi^^ag_i^^ + H.c.) where A is the superconductor 

gap and a^^^ (asko-) is the creation (annihilation) oper- 
ators in the superconducting lead with the momentum 
k = {kx,ky). Here we consider the general s-wave su- 
perconductor. The coupling Hamiltonian He is: He — 
Z^iC^^L^t' '^•S'ixi)ts*i + H.c, where the operator osi^o- = 

' tsa tsb 
-tsa ~tsb 

the parameters tsa and tsb are the coupling strengths 
between the superconductor and HgTe/CdTe QW, which 
depends on the interface's contact potential and the qual- 
ity of the coupling in the experiment. 

By using the Green's functions, the charge current 
Ine and spin current Ise from the n-th terminal of the 
HgTe/CdTc QW flowing into the device are Ine — e(/nt + 
Ini) and Ins = {h/2){In^ - Ini), whcrei^ 

1 r ( 
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Here a =;,t for a =t,4., Lt/liE) = f{E T eK), and 
fs{E) = f{E), with f{E) being the Fermi distribu- 
tion function and Vn being the voltage of the terminal- 
n. In Eq.(2), Tnm4E) = Tr{r„,G;^r^<.G^,} and 
Tnsa{E) - Tr{r„,[G'^r,G«]^^} are respectively the 
normal transmission coefficient from the terminal-n to 
the terminal-m and to the superconductor terminal, and 
T^„,msiE) = Tr{r„,G;^r^<,G| J is the AR coefiicient 
with the incident electron from the terminal-n and the 
reflected hole going to the terminal-m. The linewidth 
function Tn/siE) = i[^n/s ^ (^ri/s)^] ^^'^ Green's 
functions G''/''{E) can be calculated from G^{E) = 
iG-iE)y ^{E- H,,n - E„ - where H,e„ 

is the Hamiltonian of the scattering region as shown in 
Fig. If (dotted region). The self-energy functions SJJ(i?) 
and S^(£') due to terminals of the TI and superconduct- 
ing lead can be calculated as in Ref . [T^ - [l^ . In the fol- 
lowing numerical calculations, we choose the parameters 
from the realistic materials:^ (1) the HgTe/CdTe QW's 
parameters are A = 36A.5meVnm, B — —dSdmeVnm^ , 
C — 0, and D — —bl2meV nm? . (2) the superconduc- 
tor's parameters are the gap energy A = ImeV. The 
lattice constant a is set to Snm and the TI-S coupling 
strengths are taken tsa — tst = t. 




FIG. 1: (Color online) (a) and (b): The AR coefficient 
vs. Ef (a) and energy E (b) in the two-terminal device, (c) 



(d), and (e): 
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Ef (c), TYr^ij, vs. the energy E (d), 
and 24. vs. Ef (e) in the four-terminal device with the 
length L = 750nm. The other parameters in panels (a)-(e) 
are: the width A'' — 300nm, coupling strength t = 3meV, 
and M = —lOmeV. (f) is the schematic diagram for the 
two-terminal and four-terminal devices. 



We flrst study the two-terminal system. Due to the 
time-reversal symmetry and the C2 symmetry around 
the y axis, the AR coefficients have the properties 
T,l,^{E) = T,l,^{E) ^ T\E) and T\E) = T\-E). 

Fig. la and lb show the AR coefficient versus the 
Fermi energy Ep and energy E of the incident elec- 
tron. Now the parameter M is set —lOmeV so that 
the HgTe/CdTe QW is in the TI phase.^ Fig.la and lb 
clearly show that exhibits a plateau with the plateau 
value equal to one as long as Ep is in the bulk gap of 
TI and the energy E within the superconductor gap A. 
Since there is only one transmission channel (i.e. the edge 
state) when Ep is in the bulk gap, = 1 means that 
the incident electron is completely Andreev reflected. It 
is remarkable that the quantum Andreev effect occurs 
in such a wide range oi Ep and E, while in all previ- 
ous metal-superconductor hybrid systems the resonant 
condition T"^ = 1 occurs only at certain set of parame- 
ters. Fig.la also shows that when Ep is out of the bulk 
gap, there is no plateau in the AR coefficient and 
depends on both Ep and E. For Ep < M, can 
be larger than one (see Fig.la) because there are many 
transmission channels. However, the AR coefficient for 
each transmission channel is still much smaller than one. 

In order to reveal the nature of quantum Andreev ef- 
fect, we will study the four-terminal device. In the four- 
terminal device, the AR coefficient has 3x3 = 9 
elements, so the trajectories of the incident electron 
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FIG. 2: (Color online) T^^^ vs. t (a) and width (b) in 
the four-terminal device with Ef = —5meV, E = O.lmeV, 
M = —lOmeV, and length L — 750nm. 



and reflected hole can clearly be shown. When the 
HgTe/CdTe QW is in the TI regime, only T^^ is 
nonzero and other eight elements are zero. In particu- 
lar, as soon as the energy \E\ < A and M < Ep < —M, 
the quantized plateau with T-^ = ^ emerges and this 
plateau is independent of Ep and E (see Fig.lc and Id). 
These results can be understood with the help of helical 
edge states, in which the spin-up and spin-down carriers 
move along the edge of the HgTe/CdTe QW in clock- 
wise and counter-clockwise directions, respectively.— We 
consider the case of the spin-up electron coming from the 
terminal-1 when the energy < A. As shown in Fig. If, 
two reflection processes occur at the TI-S interface: (1). 
it can be Andreev reflected back as a spin-down hole to 
the same terminal which will contribute to j^. (2). It 
can also be normal reflected as an electron (spin up) along 
the TI-S interface, eventually to terminal-3. Note that 
the normal reflection as an electron back to the terminal- 
1 is prohibited by the time-reversal invariance and the 
helical edge states being a pair of Kramer states. The 
reflected electron traveling to terminal-3 has to go along 
the TI-S interface since the only available transmission 
channel is the edge state. This results in a reflection 
again at the TI-S interface where part of the electron is 
Andreev reflected as the hole back to terminal-1 and the 
rest is normal reflected as electron towards terminal-3. 
As this continues, multi-reflections occur as normal elec- 
tron traverses along the TI-S interface and eventually the 
transmission probability T13 becomes zero if the TI-S in- 
terface is long enough. Clearly, it is this multi-AR that 
gives rise to the quantum Andreev effect with the quan- 
tized AR T{| -j^^ — 1. Furthermore, we have three obser- 
vations: (i) If the energy \E\ > A, all AR coefficients de- 
crease as usual (see Fig. Id) because of the occurrence 
of the normal tunneling from TI to the superconductor, 
(ii) When Ep is out of the bulk gap, all AR coefficients 
have the same behavior: AR coefficient for each trans- 
mission channel is in general small and strongly depends 
on Ep and E (see Fig.lc, Id, and le). (iii) For the spin- 
down incident carrier, it is easy to show that has 
the quantized plateau similar to T-f| 

Next, we study how the quantized AR is affected by 
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FIG. 3: (Color online) (a), (b), and (c) are Ti^,^ (a), T^,2i 
(b), and rms{T^ -i^) (c) vs. Ep with the energy E — O.lmeV 
for different disorder strengths W. (d) is T^ -^^ vs. E with 
Ep — 0. The length of the disorder region Lw = 150nm and 
the other parameters are the same as Fig.l. Here all curves 
are averaged over up to 2000 random configurations. 



the system parameters. Fig. 2 shows the AR coefficient 
^it.U versus the width N of the HgTe/CdTe QW ribbon 
and the coupling strength t. The results show that the 
quantization of AR persists for a broad range of the cou- 
pling strength t. In addition, the wider the width N, the 
broader the quantization plateau is. For N — lOOOnm, 
the AR quantization plateau can sustain when t varies 
nearly one order of magnitude. Finally, with the in- 
crease of the width N, the AR coefficient -^^ rises 
monotonously before it reaches quantized value due to 
the fact that the incident electron has more chance of 
multi-AR for the longer TI-S interface. These results 
show the universal feature of the quantum Andreev ef- 
fect: it is independent of the system parameters. 

This universal feature can also be understood analyti- 
cally from multi-AR along the interface of TI and super- 
conductor. Due to the nature of TI the incoming spin 
up electron can be Andreev reflected as spin down hole 
with AR amplitude va and normal reflected as a spin up 
electron with reflection amplitude r. Note that and 
r, respectively, play the role of reflection and transmis- 
sion amplitudes for normal system. Clearly, the normal 
reflection amplitude that plays the role of transmission 
coefficient for normal system will vanish after N succes- 
sive ARs for large N. While for a particular AR ta and r 
depend on system parameters such as Ep, incoming elec- 
tron energy, and coupling strength t, it is the multi-AR 
that makes the quantum Andreev effect independent of 
system parameters. 

Is the quantized AR robust against the disorder? To 
answer this question, we consider the on-site Anderson 
disorder in a region near the TI-S interface (see the light 
gray (red) region in Fig. If). Because of the disorder, an 
extra on-site term ^I'jU'i^I'i is added on each site i in the 



4 




V (mV) V (mV) 



FIG. 4: (Color online) (a) is the conductance Gie (solid 
curves) and spin conductance Gis (dotted curves) vs. the 
bias V for different Ep (with the legend being the same in 
(b)) for the four-terminal device, (b) is the conductance Ge 
vs. the bias V for the two-terminal device, (c) and (d) are 
^it,ij. ^'^'^ conductance Gie vs. bias V (d) for the 

four-terminal device with the positive AI — 2meV. The other 
parameters in (a)-(d) are the same as Fig.l. 

disorder region, where wi is the 4x4 diagonal matrix 
with diagonal elements {wi, Wi, —Wi, —Wi). Wi is assumed 
uniformly distributed in the range [—W/2, W/2] with the 
disorder strength W. Fig. 3 shows the AR coefficients 
i4,/2^ ^^'^ fluctuation versus Ep and E. The re- 
sults show that the quantized AR plateau in -^^ is 
very robust: the AR plateau can persist and its fluctua- 
tion rms{T^ is zero for the disorder strength W up to 
lOOmeV because of the helical edge states being very ro- 
bust. Upon further increasing of W from lOOmeV, -^^ 
starts to decrease and the fluctuation becomes nonzero 
because at large disorders the system reaches the dif- 
fusive regime and the helical edge states are destroyed. 
Hence, as long as the edge state is survived, disorder has 
no effect on the quantum Andreev effect. 

Let us investigate the conductance Gne [Gne = 
dIne/dV) and spin conductance Gns {Gns = dIns/dV). 
We set the biases of the HgTe/CdTe QW terminals, 
Vi = V2 = V3 = V, and the superconductor-terminal 



bias Vs = 0. Fig. 4a shows the conductance Gie and 
spin conductance Gis versus the bias V for the four- 
terminal system. When Ep is inside the bulk gap, Gie 
(in the unit e^/h) is exactly equal to Gis (in the unit 
e/47r) since only the spin- up electron traverses from the 
terminal-I to the TTS interface with the spin-down hole 
Andreev reflected back. In particular, when V < A/e, 
a (spin) conductance plateau emerges with the plateau 
value 2e^/h {e/2n) because of the quantized AR with 
Tj^ — 1. Here we emphasize that like the quantized 
AR this conductance plateau is also universal (i.e. it 
is independent of the system parameters and the qual- 
ity of TI-S coupling) and robust against the disorder. 
Since the TI phase in HgTe/CdTe QW has been real- 
ized experimentally,—'^ this predicted quantum Andreev 
effect with quantized conductance plateau should not be 
difficult to observe using the present technology. On the 
other hand, when Ep is out of the bulk gap, both Gie 
and Gis are not equal and they are sensitive to the sys- 
tem parameters. Finally, we comment on following two 
points: (i) If HgTe/CdTe QW is in the normal state (i.e. 
M > 0), no plateau emerges for both the AR coefficient 
^nt ml ^^'^ conductance Gne (see Fig. 4c and 4d), similar 
to the ordinary conductor-superconductor hybrid system, 
(ii). For the two-terminal device with the HgTe/CdTe 
QW in TI regime, the conductance Ge also exhibits the 
quantum Andreev effect with quantized AR conductance 
Ae^/h at y < A/e (see Fig. 4b). However, the spin con- 
ductance Gs is exactly zero due to the fact that the left 
and right edges of the HgTe/CdTc QW ribbon carry the 
same current but opposite spin current. 

In summary, we predict a quantum Andreev effect in 
the 2D TI-S hybrid system, in which the AR coefficient is 
quantized with the value one. Importantly, the quantized 
AR plateau is independent of various system parameters 
and the quality of TI-S coupling, and it also is robust 
against the disorders. Hence it should be easy to observe 
it using the present technology. Due to the quantized 
AR, the conductance also exhibits the same quantized 
plateau when the bias is within the superconductor gap. 
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